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The first global quantum simulation of semiconductor-based quantum-cascade lasers is presented. 
Our three-dimensional approach allows to study in a purely microscopic way the current-voltage 
characteristics of state-of-the-art unipolar nanostructures, and therefore to answer the long-standing 
controversial question: is charge transport in quantum- cascade lasers mainly coherent or incoherent? 
Our analysis shows that: (i) Quantum corrections to the semiclassical scenario are minor; (ii) 
Inclusion of carrier- phonon and carrier-carrier scattering gives excellent agreement with experimental 
results. 
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Since the seminal paper of Esaki and Tsu [Q, 
semiconductor-based nanometric heterostructures have 
been the subject of an impressive theoretical and exper- 
imental activity, due to their high potential impact in 
both fundamental and applied research One of the 

main fields of research focuses on exploiting "band-gap 
engineering" , namely the splitting of the bulk conduction 
band into several subbands, to generate and detect elec- 
tromagnetic radiation in the infrared spectral region, as 
originally envisioned by Kazarinov and Suris 

Unipolar coherent-light sources like quantum-cascade 
lasers (QCLs) [D are complex devices, whose core is 
a multi-quantum-well (MQW) structure made up of 
repeated stages of active regions sandwiched between 
electron-injecting and collecting regions. When a proper 
bias is applied, an "electron cascade" along the subse- 
quent quantized-level energy staircase takes place. QCLs 
are usually modelled in terms of n-level systems [^. As 
pointed out in Q , such a macroscopic modeling can only 
operate as an a posteriori fitting procedure. In con- 
trast, for a detailed understanding of the basic phys- 
ical processes involved, a fully three-dimensional (3D) 
description is needed. More specifically, two main is- 
sues need to be addressed: (i) the nature of the hot- 
carrier relaxation within the device active region; (ii) 
the nature — coherent versus incoherent — of the phys- 
ical mechanisms governing charge transport through 
injector/active-region/collector interfaces. 

Point (i) has been recently addressed in , where the 
usual macroscopic treatment of the device active region 
has been compared to a fully kinetic description, based 
on a Monte Carlo (MC) solution |^ of the following set 
of equations: 



^/ki^ — [ffki/ — Tky/kiyli/i +^[-Pki/,k'i/'/kV — ^k'i/',kt//ki/] 

(1) 



phenomenological level — injection/loss (i/1) of carriers 
with parallel or in-plane wavevector k in subband v, while 
the last ones describe intra- as well as inter-snbband in- 
and out-scattering processes {k.v kV). As reported 
in 1^], the quantum-cascade within the active region is 
mainly governed by LO-phonon emission. However, such 
a microscopic analysis, being limited to the device ac- 
tive region only, does not allow to answer point (ii); This 
issue is intimately related to the long-standing controver- 
sial question ||^: is charge transport in quantum- cascade 
lasers mainly coherent or incoherent? 

To provide a definite answer to this fundamental 
question, we present the first global quantum simula- 
tion — injector plus active region plus collector — of 
semiconductor-based QCL structures. To this end, two 
basic steps are needed: First, the partially phenomeno- 
logic model in (|l|) has to be replaced by a fully micro- 
scopic description of the whole MQW core structure; Sec- 
ond, the semiclassical or Boltzmann-like treatment will 
be replaced by a fully quantum-mechanical one. These 
are not trivial tasks. 

The first step requires a proper simulation scheme to 
"close the circuit" without resorting to phenomenolog- 
ical i/1 parameters. To this aim, given the set of 3D 
single-particle electron states {ki^} corresponding to a 
single QCL stage, we consider the ideal MQW structure 
obtained as infinite repetition of this QCL periodicity 
region (see Fig. |l|). Within such extended scheme, the 
time evolution of the carrier distribution function /ka is 
governed by the following Boltzmann-like equation: 



ka.k'i 



'/k 



ka I 



(2) 
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Here, the first two terms describe — still on a partially 



Here, a = (A, v) denotes the generic electron state in our 
• MQW structure, i.e., the i^-th state of the A-th stage. 
To "close the circuit" , we impose periodic boundary con- 
ditions limiting the inter-stage (A' ^ A) scattering to 
just nearest-neighbor coupling (A' = Ail). In view of 
the translational symmetry, we are allowed to simulate 
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carrier transport over the central — i.e., A = — stage 
only 0. 

We have applied the above — still semiclassical — 
global-simulation scheme to state-of-the-art QCL struc- 
tures. As prototypical device, we have considered the 
GaAs/(Al,Ga) As-based diagonal-configuration QCL in 
PI , schematically depicted in Fig. [T], in which our sim- 
ulation strategy is also sketched. Here, the energy levels 
and probability densities of various electron states within 
the simulated stage (A = 0) are plotted: They corre- 
spond to the device active region {v = 1,2,3 according 
to the standard notation) as well as to the collector re- 
gion {i/ = A,B,C,D,E). In order to properly model phase- 
breaking hopping processes, in addition to carrier-optical 
phonon scattering, all various intra- as well as intersub- 
band carrier-carrier interaction mechanisms have been 
considered 

As a starting point, wc have investigated the relative 
weigth of the carrier-carrier and carrier-phonon compet- 
ing energy- relaxation channels. The time evolution of 
the carrier population in the various subbands as well 
as of the total current density are depicted in Fig. 2. 
Parts (a) and (b) report the population dynamics with- 
out and with two-body carrier-carrier scattering, respec- 
tively. In our "charge-conserving" scheme, we start 
the simulation assuming the total number of carriers 
to be equally distributed among the different subbands; 
then the electron distribution functions evolve accord- 
ing to Eq. (|^) and a steady-state condition is eventu- 
ally reached, leading to the desired 3 — > 2 population 
inversion. As we can see, the inclusion of intercarrier 
scattering has significant effects: It strongly increases 
inter-subband carrier redistribution, thus reducing the 
electron accumulation in the lowest energy level A and 
optimizing the coupling between active region and in- 
jector/collector (the populations of subbands 3 and B 
get equal). This effect comes out to be crucial in de- 
termining the electron flux through the MQW structure. 
Figure 2(c) shows the simulated current- voltage charac- 
teristics of our QCL device, obtained with and without 
carrier-carrier interaction. At the threshold operating 
parameters estimated in fll|| [marked by an arrow in 
Fig. 2(c)] the current density in the presence of both 
electron-phonon and electron-electron scattering mech- 
anisms is about 4 kA/cm^. This value is in relatively 
good agreement with experiments IT^ ]. 

The results plotted in fig. || clearly demonstrate that 
within a purely semiclassical picture the electron-phonon 
interaction alone is not able to efficiently couple the in- 
jector subbands to the active region ones: While carrier- 
phonon relaxation well describes the electronic quantum 
cascade within the bare active region [0, carrier-carrier 
scattering plays an essential role in determining charge 
transport through the full core region. This can be as- 
cribed to two typical features of carrier-carrier interaction 
— compared to the case of carrier-phonon — : (i) this is a 



long-range two-body interaction mechanism, which also 
couples non-overlapping single-particle states [see Fig. 1]; 
(ii) the corresponding scattering process at relatively low 
carrier density is quasielastic, thus coupling nearly reso- 
nant energy levels, like states 3 and B. 

So far, no quantum mechanical effects, like coherent 
resonant tunneling between adjacent states, have been 
considered. In order to see how such coherent phenom- 
ena can change the scenario presented so far, we have ex- 
tended our semiclassical simulation scheme in terms of a 
density- matrix formalism [ p^ . In the proposed quantum- 
transport approach the basic ingredient is the single- 
particle density matrix pij — (ajaj^, where (aj de- 
note creation (destruction) operators for a carrier in state 
i = ka [ p^ ] . Its time evolution is given by: 

i'j' 

(3) 

where hcoij = e, — ej is the energy difference between 
states i and j. Here, the first term describes the coher- 
ent evolution of the noninteracting carrier system while 
the second contribution describes energy relaxation as 
well as dephasing, in terms of the generalized in- and 
out-scattering superoperators F [jl^ . Equation (^ is the 
desired quantum-mechanical generalization of the Boltz- 
mann transport equation in (g). 

Analogous to our semiclassical simulation scheme, for 
the new quantum-transport formalism we can adopt the 
same periodic conditions to "close the circuit". More- 
over, since for the QCL design considered in-plane and 
along- z carrier dynamics are strongly decoupled [pT| , it is 
possible to adopt a factorization of the density matrix ac- 
cording to pw = Pkiy,k'i>' = Piyi/'/k(5k,k', where /il denotes 
the parallel or in-plane carrier distribution. This allows 
us to obtain an effective equation of motion for pi,,^' : the 
latter, which has again the structure of Eq. (^), involves 
effective scattering matrices, given by an in-plane average 
of the quantities r'"/°"* in (|). 

Comparing the results obtained with the proposed 
quantum-transport approach to those of the semiclassi- 
cal global simulation scheme, we find negligible quan- 
tum corrections (of a few percents) to the stationary 
current density. In the absence of carrier-carrier scat- 
tering, we get, e.g., a 2% quantum correction to the re- 
sult at threshold reported in Fig. 2. This is due to the 
extremely small value of non diagonal density-matrix ele- 
ments Pi,=^v' (compared to the diagonal ones). The phys- 
ical interpretation of such a behaviour proceeds as fol- 
lows: Non-diagonal elements in r'"/°"* tend to maintain 
a non-diagonal density matrix also in stationary condi- 
tions. On the other hand, diagonal energy-relaxation and 
dephasing processes tend to suppress non-diagonal terms 
on the sub-picosecond time-scale. Since the average de- 
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vice transit time is of the order of several picoseconds, 
the degree of coherence, i.e., the weight of non-diagonal 
density-matrix elements, in stationary conditions is very 
small. This does not mean that coherent phenomena, like 
resonant-tunneling processes, are not present. Figure 3 
presents the first simulation of the ultrafast dynamics of a 
properly tailored electron wavepacket, within the MQW 
core region of the QCL sketched in Fig. 1. The aim is 
to focus on the injector — active-region tunneling mech- 
anisms. For this reason, the system has been prepared 
at time t = to reproduce a charge-density distribution 
fully localized in the A = injector and compatible with 
resonant tunneling into the next-stage active region (see 
(0,B)-(-|-l,3) alignment in Fig. 1). As clearly shown, the 
transient dynamics is characterized by a strong interplay 
between phase-coherence and relaxation (on a picosec- 
ond time-scale); Only at much longer times it will even- 
tually reach the stationary transport solution, in which 
incoherent sequential tunneling is the dominant interwell 
mechanism |p^ . 

In summary, we are in the position to answer the long- 
standing controversial question on the nature — coherent 
versus incoherent — of charge transport in QCLs previ- 
ously mentioned: For the typical structures considered, 
energy-relaxation and dephasing processes are so strong 
to destroy any phase-coherence effect on a sub-picosecond 
time-scale; as a result, the usual semiclassical or incoher- 
ent description of stationary charge transport is found to 
be in excellent agreement with experiments. 
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FIG. 1. Schematic representation of the conduc- 

tion-band profile along the growth direction for the diago- 
nal-configuration QCL structure of Ref. The MQW is 
biased by an electric field of 48 kV/cm. The levels = 1, 2, 3 
and 1/ — A,B,C,D,E in the active and collector regions (full 
and dashed lines, respectively) of the simulated stage (A — 0) 
are also plotted together with the corresponding probability 
densities. The replica of level 3 in the following stage A = -1-1 
is shown for clarity. 
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FIG. 2. Time evolution of simulated carrier densities 
in the various subbands of the MQW structure of Fig. |l| 
(full lines: v = A,B,C,D,E; dotted lines: v = 1,2,3), with- 
out (a) and with (b) intercarrier scattering, (c): simulated 
applied-field vs current-density characteristics of the whole 
structure at 77 K, in presence (discs) and absence (filled 
squares) of carrier-carrier interaction. Threshold applied field 
(48 kV/cm) is marked by an arrow. Dashed lines are a guide 
to the eye. 

FIG. 3. Time evolution of the charge density for an elec- 
tron wavepacket properly tailored to study the carrier tun- 
neling dynamics accross the injection barrier, for the QCL 
design of Fig. 1. At t = ps (lower panel (c)) the wavepacket 
is fully localized in the injector. Shaded regions correspond to 
(Al,Ga)As barriers in the heterostructure design. The tran- 
sient dynamics (middle panel (b)) is characterized by a strong 
interplay between phase-coherence and relaxation processes. 
At much longer times (upper panel (a)) the system will even- 
tually evolve into the stationary transport solution. 
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